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ABSTRACT 

Data obtained by the Skylab radar altimeter (S-193) 

over ”the' _ ContirnentaT USA ‘ for Pass T of'S Ii^2 — and “P aTss e s TF7 

17 j_ 18* and 28 of SL-3 were analyzed by the Naval Research 
Laboratory to relate radar return signal parameters such as 
range, power, and waveform to ground track terrain features. 
The altimeter, which was originally designed for ocean 
observations, performed predictably well over terrain with 
moderate to low relief. The maximum radar return from more 
complex areas corresponded usually to the lower flatter 
regions within the 7-kilometer diameter footprint. 

The altimeter-measured heights correlated with map 
topography with an rms deviation of ±30 meters for terrain 
topography with large height dispersion (Passes 1 and 28) 
and ±5 meters for slowly varying terrain (Passes 16, 17, and 
18) . The received power varied with specific terrain 
features and it ranged from a radar cross-section of 100 dB- 
above-one- square-meter (dBsm) for the Utah Salt Flats to 
60 dBsm for mountains with forests. The ocean radar cross- 
section was located approximately at the mid-point of these 
extremes (80 dBsm) . 

The analysis of the waveform and statistical character- 
istics of the individual signals indicated that larger 
signal power was usually related to dominantly specular 
reflections from patches 100 meters or larger. 

The results established the capability of a spacecraft 
altimeter system to observe characteristic radar signatures 
for different types of terrain within the bandwidth restric- 
tions of the AGC and range-tracker loops. Improved reliabi- 
lity of operation and more complete height information on 
complex regions can be attained by a dedicated altimeter 
system with extended dynamic range in both power and altitude* 
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SKYLAB ALTIMETER OBSERVATIONS OVER TERRAIN 


INTRODUCTION 

The narrow-pulse radar altimeter of the Skylab EREP 
package (S-193) was operated over land areas to determine 
from satellite heights the radar reflection characteristics 
of the terrain at nadir, and to evaluate the ability of a 
satellite radar altimeter to sense the topographic and 
physical properties of terrain reflecting surfaces. Obser- 
vations were made over different types of terrain of the 
continental USA and the range, received power, vmveform, and 
statistical properties of the radar return signal as a 
function of the sub-satellite area were analyzed. Although 
the skylab altimeter was originally designed to operate over 
the ocean surface where relatively homogeneous reflection 
conditions prevail, it performed predictably well over 
terrain with moderate to low relief,* with range-tracker loss 
only over land areas with large variations and complex 
reflection properties. 

In this report are described pertinent technical 
parameters of the Skylab (S-193) radar altimeter system and 
its operational modes, which determine the types and the 
character of the data obtained (Section 1) , the data pro- 
cessing and analysis undertaken at the Naval Research 
Laboratory (Section 2) , error analysis (Section 3) , the 
results of and the conclusions drawn from the data analysis 
(Sections 4 and 5) , and consideration of parameters for a 
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dedicated radar altimeter system with improved response and 
data handling capability (Section 5) . 

1. ALTIMETER DESIGN AND OPERATION 

In this section are given the basic engineering param- 
eters of the Skylab (S-193) radar altimeter system, the 
general operational geometry of radar altimetry, the details 
of operation of the automatic gain control (AGC) , the range 
acquisition, the range tracker, and the sample— and— hold 
(S&H) gate wave-sampling subsystems of the altimeter, and 
also information concerning the pre-programmed format used 
for data acquisition and calibration. Emphasized are the 
design and operational capabilities and some of the limita 
tions of the altimeter to determine land relief and asso- 
ciated terrain features suitable for cartographic applications 

1.1 Geometry of Terrain Observations 

The basic Skylab altimeter system parameters are given 
in Table 1. These parameters determine the operating 
geometry shown in Figures 1 and 2. Two basically different 
types of terrain can be considered, each with its character- 
istic radar signal return. In Figure 1, a flat area is 

assumed and the resulting radar return may consist of either 

% 

a specular or diffuse reflection or a combination of both. 

If the diffuse reflection is dominant, the return will be 
similar to that obtained from the ocean, as shown in the 
lower right. In this case, the mean radar return consists 
of a rise-time of about 100 ns, corresponding to the 100-ns 
transmitted pulse, and a trailing edge of about 500 ns 
produced by the antenna drop-off as the incident angle 
increases. The radar illuminated footprint in the first 
100 ns is of the order of 7 km and will increase to about 
10 km in the trailing edge. For specular reflection, the 
radar return waveform is approximately the same as that of 
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TABLE 1 


NOMINAL SKYLAB ALTIMETER SYSTEM PARAMETERS 
FOR TERRAIN OBSERVATIONS 


Transmitter 

TWT Power Output 
Frequency 
Pulse Width 

Pulse Repetition Frequency 


2 kW 
13.9 GHz 
100 ns 
250 pps 


Antenna 

Paraboloid Diameter 
Aperture Efficiency 
Gain (X = 2.16 cm) 

3 dB beamwidth 


1.12 m 

50% 

41.3 dB 
1.6° (measured) 


Receiver 

Pre-Amplifier (Tunnel Diode) 
IF Center Frequency 
Bandwidth 


33 dB 
350 MHz 

10 and 100 MHz 


Derived Radar Footprint Diameter 
for Skylab Height (4.3 x 10 m) 

Specular Return 
Pulse Width 
Beamwidth (3 dB) 


100 m 
7.3 km 
8 . 5 km 
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the transmitted pulse as shown in the upper right of Figure 2. 
The reflecting area in this case may correspond to the 
maximum area that allows a phase coherent return, and 
which, at the satellite height and the RF wavelength used by 
the radar, has a diameter of the order of 100 m. In general, 
the radar return will be composed of both diffuse and 
specular reflections. If more than one smooth area within 
the radar footprint contributes to specular reflection, a 
fluctuating signal with a pulse shape similar to that of the 
transmitted pulse is produced. 

In mountainous regions where the vertical height struc- 
ture of the observed terrain may extend over several hundred 
meters within the footprint, the waveform of the radar 
return may appear as shown in Figure 2. Specular reflection 
may be present at different heights in a rough area such 
that return delay times from such reflections will extend 
over several microseconds and will also vary due to different 
scattering characteristics. It will be shown later that the 
AGC operation of the altimeter system will favor the larger 
signals for range lock operation, and that only a small 
portion (250 ns) of the return is utilized for the radar 
Output measurements. In mountainous regions the spot size 
will be limited by the antenna beamwidth rather than by the 
pulse width and it will thus be of the order of 10 km. 

1.2 Basic Skylab Altimeter Operation 

The Skylab radar altimeter is designed to measure the 
amplitude and range of a received signal with high precision. 
The basic signal operations performed by the receiver 
system, shown in Figure 3, include the following functions: 

1. Automatic gain control (AGC) of the receiver, 

2. Range acquisition, 
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Determination and tracking of range by the range 
tracker, 

4. Reproduction of a portion of the received waveform 
by eight 25-ns sample-and-hold (S&H) gates. 

The dynamic range of the detector is about 6 dB, and 
since the dynamic range oi the received signal is expected 
to be considerably larger, the AGC establishes the proper 
power level at the detector so that the other three functrons 
can be performed with a signal level that is withrn the 
linear range of the detector. The acquisition process 
narrows the range uncertainty so that the range tracker is 
able to lock and track within a small time delay region of 
the radar return. Once range lock is verified, the AGC gate 
is matched to the expected received waveform, the S&H gates 
are positioned in a fixed time sequence relative to the 
range reference output of the range tracker, and the system 
initiates the pre-programmed sequence of data acquisition. 

The data outputs of the altimeter are recorded at fixed time 
intervals and include the AGC (received power) , the time 
delay of the transmitted radar pulse (range) , and the eight 
S&H gated amplitudes (waveform and statistical amplitude 
fluctuations) . 

The signal processing characteristics are summarized m 
Table 2 and will be considered in more detail in the follow- 
ing sections . 

1.3 AGC Operations 

The automatic-gain-control (AGC) circuit senses the 
power level of the radar return signal within a pre-positioned 
time gate in which the signal is expected, and it then 
utilizes the mean of 64 peak signal amplitudes to control 
the gain of the receiver. Before acquisition an 800-micro- 
second wide sampling gate and following acquisition and 
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TABLE 2 


SIGNAL PROCESSING SYSTEM PARAMETERS 


AGC 


Acquisition Window 
Tracking Window 
Closed Loop Time Constant 
Dynamic Range 
.AGC Output 
AGC Granularity 
Steady State Tracking Error 


800 ys 
600 ns 
.2 sec 
45 dB 

Peak of 64 pulse average 
.1 dB 
1/2 dB 


Altitude Tracking Loop 


Type 

Loop Bandwidth 
Altitude Output 
Altitude Granularity 


Digital, 200 MHz logic 
2 Hz 

32 pulse average 
2.5 ns 


Waveform Sampling 


Number of S&H Gates 
Sampling Gate Width 
Gate Spacing 


8 

25 ns 
25 ns 


during tracking a 600-nanosecond wide gate positioned at the 
tracked range are used to gate the return signal for the 
AGC. 

To allow for the normalized cross-section a Q of differ- 
ent sea states, the height variations of the spacecraft, 
attitude changes, and different altimeter mode operations, a 
large dynamic range for the AGC is provided as shown in 
Figure 4, which extends from about -10 dBm to -60 dBm. The 
noise threshold, which is also shown, is approximately 
-50 dBm for the 100 MHz receiver bandwidth sub-mode and 
-60 dBm for the 10 MHz sub-mode. 

The large dynamic range of about 50 dB permits the 
altimeter to operate over a variety of terrain. However, 
Since the relatively rapid changes in received power level 
encountered in observations over solid terrain are attenuated 
by the AGC smoothing filter, which was designed for the 
gradual changes in ocean reflections, rapid changes in 
received power will cause loss of range lock, as well as 
reduce the acquisition capability of the altimeter. 

1,4 Range Acquisition 

Range acquisition is achieved by performing two thres- 
hold decisions respectively in (1) the coarse-search mode 
and (2) the fine-search mode. A graphic description of the 
two search modes is given in Figure 5. In the coarse- 
search, a one-microsecond gate scans in 320-ns steps at a 
rate of 4 ms/step for an interval of 82 microseconds . 

(12,300 m) which is centered about an a priori estimated 
time delay T . The received signal from the sub-satellite 
area, which is shown to the left of the estimated arrival 
time, is nominally a non-symmetrical pulse with about a 
100-ns rise-time and a 500-ns trailing edge. Since the 
maximum rate of change in the delay time of the received 
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signal, due to orbital vertical motion and ocean surface 
deviations from a sphere, is of the order of a few nanoseconds 
per 4 ms, the received signal is essentially frozen relative 
to the 320-ns shifts of the one-microsecond gate. Thus, as 
the -search gate scans across the signal, there are at least 
2 to 3 gates that will contain the received signal. When 
the signal is detected more than once, as indicated by a 
threshold decision, the acquisition process switches to the 
fine-search mode. This mode uses a 200-ns gate in 50-ns 
shifts at the 4-ms/step rate, which scans a 3 . 2 -microsecond 
time interval that is centered at the detected coarse-search 
time delay. The same decision threshold process is repeated 
and when the signal is located more than once in the fine- 
search mode, the range tracker operation is initiated. If 
the signal is missed in the acquisition process, the opera 
tion reverts back to the coarse-search mode and the scan is 
repeated. The maximum time of acquisition, assuming no miss, 
is 1.04 seconds in the coarse-search mode, 0.256 seconds m 
the fine-search mode, and about 0.5 seconds for range 
verification for a total of 1.8 seconds. 

The acquisition process over the ocean is based on the 
assumption that the received signal has a unique and narrow 
pulse shape, constant mean power, and is located at an 
approximately fixed time delay inside the a prion time 
delay uncertainty. The acquisition process over terrain may 
be considerably altered because: 

1. the pulse shape will, in general, be extended in 
time and vary with different types of reflection 
characteristics and terrain topography, 

2. the mean power may change relatively rapidly at 
different portions of the return, and 
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3. the location of the ramp portion of the return, on 
which the acquisition system operates, may shift 
rapidly during the time of detection. 


A possible radar return from a mountainous area is 
shown in Figure 2. Acquisition here is possible at several 
points where the signal amplitude increases over a few 
hundred nanoseconds. However, the earlier larger signals 
are favored because the AGC is controlled by the peak 
amplitude of the return and the scan is from left to right. 
Thus, several possibilities exist which could prevent 
efficient range acquisition. For example: 

1. The coarse-search may detect a signal at an early 

• point in time that does not coincide with the peak 
signal detected by the AGC. The fine-search may 
not be able to detect the signal at the same 
location with AGC loop control overly reducing 
appropriate receiver gain, or on the other hand, 
the fine- search can detect the signal but during 
pull-in the range tracker is not able to lock in 
because of the small available signal. In either 
case, the acquisition reverts to the coarse-search 

mode and a new attempt is made to locate a signal 

% 

ramp . 

2. The coarse-search may detect a signal that is the 
peak signal, but the mean amplitude of the signal 
may decrease rapidly due to variations in the 
reflection mechanism. The slowly responding AGC 
is not able to follow the drop in the signal level 
and neither the fine-search nor range tracker will 
be able to locate the signal. 

The coarse-search detects a signal as in paragraph 
2 above, but the signal moves rapidly in range in 
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a direction opposite to that of the scanning 
operation. Again either the fine-search will miss 
the signal location or the range tracker will not 
be able to follow the rapid range change. 

Any combination of the above conditions may increase 
the acquisition time and thus, it is not surprising that in 
areas of complex relief, acquisition times of the order of 
10 seconds or more actually occur. For some extremely 
complex terrain, the signal is not acquired at all. 

1.5 Range Tracker Operation 

Once range acquisition is achieved, the three main 
operations of the receiver system, namely, AGO, range 
tracking, and S&H gate sampling, are performed on a few 
hundred nanosecond segments of the acquired signal. Thus, 
for complex topographic terrain observations, the receiver 
system ignores any radar return whose time delay falls 
outside the narrow time interval for range tracker process- 
ing. The range tracker operation is controlled by three 
100 -nanosecond gates, as shown in Figure 6. The noise gate 
is placed 150 ns before the expected ramp of the radar 
return, and the plateau gate is placed 50-ns after the ramp 
gate. The range is tracked in a feedback loop that utilizes 
in a split gate configuration the difference between the 
amplitude of the ramp and 1/2 the amplitude of the plateau 
gate to establish the position of the radar return relative 
to the position of the split gates as shown in Figure 7. 

This voltage difference is applied to the time discriminator 
control curve m the feedback loop, which shifts the split 
gates to a position that minimizes the output voltage. 

Since the radar return is both pulse and beamwidth limited 
as indicated by the return pulse wave shape in Figure 7, the 
null of the time discriminator curve does not occur at the 
origin. 
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For terrain observations, range tracker operation may 
be modified and complicated by several factors. First, loss 
of range verification may be produced by energy in the noise 
gate reflected from higher elevations. Second, the waveform 
of the radar return, because of the different reflection 
mechanisms that may be encountered, may be considerably 
different from that assumed for ocean returns. This will 
introduce both an error in the range measurements as well as 
affect the range tracking capability. Finally, rapid 
changes in amplitude or range of the sampled radar return 
that cannot be followed by the AGC or range tracker loops 
respectively will cause either marginal operation or loss of 
range lock. 

1.6 S&H Gate Operation 

Portions of the radar return waveform used in the split 
gate configuration are sampled by eight 25-ns contiguous 
gates that are shifted at about 20-second intervals across 
the nominal noise, ramp, and plateau regions to cover 600-ns 
of the waveform shape, as shown in Figure 6. The gate 
outputs, recorded at a rate of 100 amplitudes per second, 
can be used to study the statistical characteristics of the 
radar return as a function of time delay, assuming a range 
jitter of less than 25 ns. 

Over terrain, it is not possible, in general, to obtain 
the complete waveform of the received signal, since the 
waveform may extend over several microseconds. In addition, 
the noise, ramp and plateau portions of the waveform cannot 
be combined, since the terrain and reflection character- 
istics of the observed areas may change considerably during 
the one-minute interval required to sample the 600-ns 
portion of the received waveform. However, an estimate of a 
small segment of the waveform can be made over uniform 
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terrain by combining the noise-and-ramp, or ramp-and-plateau 
portions of the waveform when the position of the sampling 
gates changes. This changeover occurs three times in each 
run, as will be shown in the next section. 

1.7 Pre-Programmed Sequence of Data Acquisition 

The altimeter measurements are grouped in a sequence of 
frames, with each frame containing data for a 1.04-second 
interval, corresponding to 7.5 kilometers of ground track 
distance or covering an extended footprint area roughly 7.5 
x 15 kilometers. The sub-satellite area is probed with a 
100-ns pulse 256 times during each frame and the radar 
system processes the received signals to produce 4 AGC 
outputs, 8 range outputs and 104 amplitude outputs for each 
of the 8 S&H gates. 

The programmed sequence for mode 1, which was exclu- 
sively used for terrain observations, is shown in Table 3. 
The mode is divided into six sub-modes (SM) . The first 
three sub— modes are used for data acquisition with the first 
sub-mode operating with a 10— MHz receiver bandwidth, the 
second sub— mode using a 100— MHz bandwidth and the antenna 
being tilted by about 0.5° in pitch angle in the final sub- 
mode. The length of each sub-mode in terms of number of 
frames and actual time duration is also given. 

In the second set of three sub-modes, the conditions of 
the first two sub-modes are simulated, and the transmitted 
pulse shape, transmitted power and internal system delay are 
calibrated. The total run time for mode 1 is about 3.5 
minutes. - — - - ~~ — — 

Each of the data acquisition sub-modes is further 
divided into three sub-sub-modes (SSM) to shift the time 
position of the S&H gates across the noise, ramp, and 
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plateau regions of the received waveform. These shifts 
occur, as indicated in Table 3, at about 15 or 20 second 
intervals, which correspond to 112 and 149 kilometers of 
ground track distance respectively. Geographic areas 
covered in the different sub-modes are listed in Tables A-l 
through A- 5 in the Appendix. 

For runs over complex terrain where frequent range 
tracker loss occurred not all sub-modes were activated. In 
general, sub-modes 1 and 2 were completed but the correct 
number of available data points was not obtained, due to 
improper operation of the frame counter. 
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*15- or 20-frame series represent three sub-sub-modes (SSM) , 
where S&H gates are positioned sequentially at noise, ramp, 
and plateau sections of signal. 


2. DATA FLOW, REDUCTION, AND ANALYSIS 

In this section are outlined the overall organization 
of Skylab altimeter data flow, and the methods of data 
reduction and analyses applied to these data at NRL to 
evaluate altimeter performance and capabilities. The 
analysis is basically a comparison of reduced altimeter data 
and ground truth, the latter mainly pertinent topographic 
information from 1:250,000 scale maps (AMS edition), supple- 
mented by some Skylab photography and weather data. Appro- 
priate references are made to detailed descriptions and 
discussions of typical cases appearing later in Sections 4 
and 5. 

* * 

2.1 Overall Data Flow 

The overall data flow from Skylab altimeter output to 
NRL is shown in Figure 8. The altimeter output was recorded 
in the spacecraft in Pulse Code Modulation (PCM) counts on 
raw data tape. The format of the altimeter output on the 
raw data tape was prescribed for a data frame, covering 1.04 
seconds in spacecraft time, 7.5 kilometers in ground track 
distance, or a footprint area roughly 7.5 by 15 kilometers, 
with 1,040 ten-bit words read out at a rate of 10 kilobits/ 
sec, giving frame time, four AGC values, eight range delay 
times, 104 sets of eight S&H gate outputs, averaged noise, 
ramp, and plateau gate outputs and associated housekeeping 
data. Details of this format are available in reference 1. 

The raw data tape was brought back after each Skylab 
mission to the Johnson Space Center (JSC) where all data, 
except the S&H gate outputs, were converted to engineering 
units. In addition, the Airlock Module Time (AMT) was 
converted to GMT, and the spacecraft position and height 
were computed from the Skybet Orbit Program. The modified 
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data stream was recorded on pre-processed tape in computer 
compatible format and transmitted to NRL with a header 
record which contained the required algorithms for read out 
of altimeter data. The data format for the pre-processed 
tape is given in references 2 and 3. 

The altimeter data on the pre-processed tape was 
further processed at NRL, where a geodyne orbit was used to 
refine spacecraft position and to convert altimeter range 
measurements to ground heights. The results were then 
compared with Skylab photography and ground truth data. 

2.2 Data Reduction 

The scheme of data reduction and analysis performed at 
NRL on the pre-processed tape is shown in Figure 9. The 
data reduction consisted of a primary evaluation of the data 
on the pre-processed tape, followed by the application of 
corrections and time tags to AGC and ALT data, and finally 
the conversion of ALT to ground height. 

2.2.1 Primary Evaluation 

The following data were read out directly for each 
frame from the pre-processed tape for primary evaluation of 
data and cross indexing of frame number and spacecraft 
location: 

1. Frame or data block number 

2. Operating mode, sub-mode (SM) and sub-sub— mode 
(SSM) 

3. Frame counter 

4. Altimeter lock status 

5. Average altitude for each frame 

6. Standard deviation of altitude within frame 
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Residual height relative to Skybet height estimate 
8. Average AGC for given frame 

9/ 10,11. Average noise, ramp, and plateau gate outputs 
for given frame 

12. Frame time (GMT) 

13, 14. Longitude and latitude of sub-satellite 

position from Skybet orbit. 

The first criterion to be met was a check of range 
tracker lock. Because of the complex reflection character- 
istics of the terrain, three separate conditions had to be 
satisfied simultaneously before range lock was considered 
acceptable . 

1. The altimeter lock status word was present (digit 
0) 

2. The ratio of ramp-to-noise gate output was larger 
than two. 

3. The difference between adjacent range readings did 
not exceed 500 m. 

In addition, estimates of the quality of data could be 
obtained from the AGC values and changes in the standard 
deviation of the altitude readings. All data satisfying 
these range lock conditions were then subjected to further 
corrections, analysis, and evaluation. 

2.2.2 Altitude Corrections 

Each altitude value was corrected for the 

following : 

1. A 40-ns delay introduced by JSC during conversion 
to engineering units 
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The instrument system delay 

3. A decrease in bandwidth delay for the 100-MHz sub- 
mode 

4. The propagation delay in the atmosphere. 

For SL-3 and SL-4 data, the corrections for the 40-ns 
delay and the instrument system delay were already made by 
JSC on the pre-processed tape, and only the two remaining 
corrections had to be applied. The instrument system and 
bandwidth delay times were derived from pre-flight and in- 
flight calibration. 

The instrument delay was derived in the following 
manner : 

In-flight calibration system delay 
(Mode 1, SM5 (CDS 3), 10 MHz bandwidth) 

Pre-flight calibration (digital delay 
generator offset, 10 MHz bandwidth) 

(reference 4) 

System delay in calibration mode (10 MHz) 

Antenna hardware delay (reference 4) 

Correction for Data Acquisition' (triangular 
Pulse) minus calibration mode (square pilse) 

Total system delay (10 MHz bandwidth) 

Pre-flight calibration correction for 
100 MHz bandwidth (SM 2) (reference 4) 

Total system delay (100 MHz) 90 ns 

For SL-2, the 40-ns delay introduced by JSC was added 
to produce a total delay correction of 220 ns for the 10 MHz 
sub-mode . 


2,624,085 ns 

2,624,000 ns 
+85 ns 
. 10 ns 

85 ns 
180 ns 

-90 ns 
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Finally, an atmospheric correction of 17 ns (Reference 
8) was applied to account for the mean delay of the signal 
through the troposphere and ionosphere. Thus the basic 
range corrections applied to the measured value t are given 
by 


t = t - t . , - t. . 
c m id bias 

where 


t atm + t(10 ° MHz) 


and 


'id 


180 ns 


tbias = 40 ns 

‘ata = 17 ns 


(system delay at 10 MHz) 
(bias introduced by JSC) 
(atmospheric correction) 


t(100 MHz) = 90 ns. 

The net effective corrections applied to the measured 
time delay in SL-2, 3, and 4 are summarized in Table 4. 


TABLE 4 

EFFECTIVE TIME DELAY CORRECTIONS AT NRL 

10 MHz 100 MHz 

* 

-237 ns -147 ns 

-17 ns 73 ns 

2.2.3 Time Tag of Output Data 

A unique time tag for each altitude and AGC 
measurement cannot be defined, . since each altitude and AGC . 
value depends on the non-uniform weighting of past values by 
the range tracker and the AGC loops, and also on the input 
characteristics of the signal. It is possible however, to 
estimate the time of each measurement from the known time 


SL-2 
SL- 3,4 
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characteristics; of the range tracker and AGC processing 
loops. If a uniform input signal and uniform weighting are 
assumed, a unique measurement time, which on the average 

will correspond to the true measurement time of each output 
value, is obtained. 

Range Time Tags 

The nominal measurement times of the eight altitude 

readings in each frame relative to frame time are listed in 
Table 5 (reference 5). 


TABLE 5 


TIME TAG FOR RANGE MEASUREMENTS RELATIVE TO FRAME TIME (FT) 


Tline Readout Measurement Time 
Posxtion Relative to FT Relative to FT 
m Frame in ms in me 


1 

66 

-80 

2 

192.5 

66 

3 

320.5 

192.5 

4 

448. 5 

320.5 

5 

576 

448.5 

6 

704,5 

576 

7 

832.5 

704.5 

8 

960 

832.5 


Time 

Separation 
in ms 

Averaging 
Time for 
each 
altitude 
in ms 

146 

128 

126.5 

125 

128 

128 

128 

128 

127.5 

128 

128.5 

127 

128 

129 

127.5 

128 


The table also gives the time of readout for each altitude, 
the time separation between sampled altitudes, and the time 
interval over which each altitude output is averaged. it is 
seen that the altitudes are not read uniformly and that 
there is an interval of 19 ms near the beginning of each 
frame where altitude is not measured. The interval in which 


32 range-tracker-processed altitude measurements are 
averaged is about 128 ms in time or 919 m in ground track 

distance. . 

Based on the time distribution shown in Table 5, the 
time tag for the averaged altitude for each frame is given 

by 


fc av alt/frame 


= FT + 383 ms . 


AGO Time Tag 

Similar time tags can be assigned to the AGC measure- 
ments. Using the AGC loop time-constant of 200 ms, the AGC 
measurement times are advanced by 200 ms as shown xn Table 6 


TABLE 6 
AGC TIME TAG 


AGC Value 
Position 
in Frame 

Time of Readout 
Relative to FT 
in ms 

Measurement Time 
Relative to FT 
in ms 

Time 

Separation 
in ms 

1 

61 

-139 

320 

2 

381 

181 

280 

3 

661 

461 

240 

4 

901 

701 

200 


The measurement time for the averaged AGC value for one frame 
is 


t a , v AGC/frame 


= FT + 301 ms. 
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S&H Gate Time Tag 

The measurement times for the S&H gat e_ values have been 
— assumed to - correspond to"the~ times - of - readout~wTiich occur - at - 
10-ms. intervals starting with the FT. . 

The measurement times of all of these values are 
decreased by about 1.6 ms to account for the one-way time 
delay of the radar signal. 

2.2.4 Conversion of Altimeter Altitudes to Ground Heights 

The measured altimeter altitudes were converted to 
ground heights relative to the geoid along the Skylab ground 
track as shown in Figure 10. 

The ground height values were obtained from 
h t - (h c + Ah - h g > - h ffl 

where h is the computed spacecraft height relative to the 
c 

SAO reference surface used in the Geodyne orbit. 

Ah • is a correction for shifting the reference surface 
(SAO) used in Geodyne orbit computation to the 
reference surface of the Marsh-Vincent Geoid (1973) 

h is the geodetic deviation obtained from the 
Marsh-Vincent Geoid (1973) , 

and 

h is the measured altitude. 

m 

Since the same flattening is used in both reference 
surfaces, the value of Ah is a constant and equal to the 
difference between the SAO equatorial radius of 6,378,155 m 
and the Marsh-Vincent equatorial radius of 6,378,142 m, or 
13 m. 
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2.3 Data Analysis 

Analyses of the reduced data performed at NRL to 
evaluate Skylab altimeter performance consisted of the 
following: 

a. Comparison of fine-structure AGC output and fine- 
structure ground heights, 

b. Comparison of altimetry ground height and map 
topography, 

c. Comparison of AGC signal level and height 
residual, 

d. Sample-and-hold (S&H) gate output calibration and 
analysis, and 

e. Correlation of AGC, radar altimetry, and topo- 
graphic map profiles. 

Fine-structure data, consisting of four AGC and eight 
range-delay times per data block, were used for comparative 
analysis only in the first and last cases. In the last 
case, (e) , running averages of fine- structure data, including 
map elevation readings every one-fourth of a data block, 
were compared. 

% 

In the second and third types of analyses conducted, 

AGC and altimetry ground height data, averaged per data 
block, were compared with ground-truth, that is, corre- 
sponding terrain elevations and other features extracted 
from 1; 250, 000 scale (AMS edition) topographic maps of the 
U.S. The data frame or data block (DB) , equivalent to 1.04 
seconds of real Skylab transit time or 7.5 kilometers in 
ground track distance, was a convenient unit for production 
processing of the data stream, and also for averaging out 
large fluctuations in range and AGC data, whether due to the 
real terrain or due to tracking and AGC loop error . It was 
also noted that some average terrain feature could be 
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identified with a data block. The details of topographic 
map use are described in Section 2.3.2. 

The sample-and-hold gate output calibration and 
analysis was performed to statistically evaluate target 
reflection characteristics, namely the type of reflection, 
specular or diffuse, and the size of reflecting patches. 

All of the reduced altimeter data averaged per data 
block, topographic map data per data block, and pertinent 
geographical descriptions are collected into a set of 
tables, "Catalogue of Terrain Radar Signatures," appendixed 
to this report. 

2.3.1 ine-Structure AGC and Ground Height 

Fine-structure AGC signal levels and the reduced 
ground heights as a function of measurement time were 
plotted for selected portions of each pass to obtain an 
understanding of the interaction between the range tracker 
and AGC loops resulting from changing topography and reflec- 
tion mechanism of the terrain. Skylab photography, available 
for Pass 28, was also included to aid in the interpretation 
of terrain influence on altimeter function. The data are 
discussed in Section 4 . 2 and illustrated in Figures 13 
through 20. 

2.3.2 Comparison of Altimetry Ground Heights and Map 

Topography 

As previously mentioned, the fine-structure radar 
altimetry ground heights' were averaged over a data block for 
the purposes of this comparison. This averaging was applied 
to all good ground height data to generate a large data base 
for the detailed evaluation of overall altimeter performance. 
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With respect to identifying corresponding map topography 
different approaches in the selection and the use of map 
elevations within a footprint were used. An initial attempt 
to generate detailed two-dimensional profiles of the ground 
track by combining a series of profiles of the terrain along 
lines parallel to the ground track was found to be laborious 
and unproductive because of large errors in the altimetry 
and map reading. Another attempt considered the identifi- 
cation of a "dominant" area within a footprint, but this 
failed because the largest areas within a footprint of a 
uniform terrain type were not ‘ necessarily the most reflecting 
This procedure was used only for Passes 1 and 28, where the 
elevation of the dominant area is identified as "expected 
height from map reading" in Figures 21 and 25. 

The procedure finally adopted for general evaluation 
of the relation between altimetry and map elevations by 
data block was to determine maximum-minimum elevation bounds 
within a footprint from the topographic maps as shown in 
Figure 10, and then observing whether the altimetry ground 
heights fell within these bounds or not. This relation is 
discussed in Section 4.3 and illustrated in Figures 21 
through 25. The maximum-minimum bounds for each data block 
are given by calculating (h mid ± Ah/2) from the data tabu- 
lated in the appendixed Catalogue of Terrain Radar Signatures 

2.3.3 AGO and Height Residuals 

The deviations of the derived altimetry ground heights 
from the map height bounds derived previously are tabulated 
and also plotted separately as height residuals and compared 
with the AGC outputs. This plot relates apparent height 
errors to received power variations and it permits the 
evaluation of the instrument reponse of the range tracker 
to changes in AGC. The data are discussed in Section 4.4 
and illustrated in Figures 26 through 30. 


2 . 3.4 


Sample-and-Hold Gate Output Calibrations and 
Analysis 

The S&H gate outputs were recorded in PCM counts 
on pre-processed tape with two bits used to signify proper 
acquisition (reference 6). These values were converted to 
voltage by using calibration curves for receiver temperature 
for each gate given in reference 7. The S&H gate calibration 
curves were actually averaged for the temperature range 
-20°C to +25°C to obtain a response for the relatively 
constant receiver temperature of 0°C. In addition, the 
amplitudes of each gate were corrected for AGC changes 
within a frame with the maximum AGC value normalized to 
unity. 

The corrected S&H gate outputs were analyzed in data 
blocks of one frame each and the following operations were 
performed. 

1. Plot of selected individual waveforms. 

2. Computation and plot of the average value and 
standard deviation of each gate output for 
selected frames to obtain the pulse shape of the 
radar return. 

% 

3 . Calculation of the amplitude distribution for each 
gate within a frame to statistically evaluate the 
amplitude fluctuations. 

4 . Computation of the amplitude auto-correlation for 
each gate for shifts of several pulse periods to 
determine the coherence of the radar return signal 
in each gate. 

This analysis is discussed in Section 4.5 and illu- 
strated by Figures 31 through 41. 
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Selected information concerning the ramp or rise time 
part of the return signal can be found in the appendixed 
Catalogue of Terrain Radar Signatures. 

2.3.5 Correlation of AGC , Radar Altimetry, and topographic 

Map Profiles 

The cross-correlation between sets of smoothed 
altimetry and map elevation data was calculated to obtain a 
quantitative evaluation of altimeter performance in measuring 
terrain elevations. Running averages of the data were used 
not only to smooth out the data, but also to match resolu- 
tion scales for meaningful comparisons to be made. Details 
of these calculations and results are described in Section 
5.1 and illustrated in Figures 44 through 47. 

A correlation of the fine topography with AGC signal 
level was also attempted, but the results were inconclusive 
based on a simple terrain reflection model. However differ- 
ent mean received power levels could be associated with 
specific types of terrain. 


3. 


IDENTIFICATION OF ERROR SOURCES 


There are a number of systematic and random errors 
affecting the accuracy of the final altimeter-derived ground 
height, ranging from various spacecraft position and time 
errors, which are summarized in Table 7, to errors in the 
measurement of range and AGC signal power, and the sampling 
of radar return signal waveform. The measurement errors are 
of the order of 100 to 200 meters in range, and they are 
mainly due to poorer tracking performance over terrain with 
high complex relief. 

Errors in reading topographic map elevations were of 
the order of 10 to 30 meters, assuming that the correct 
corresponding ground track was known. 

3.1 Position Uncertainty of Altimeter Measurements 

The accuracy of the location of the reflecting area 
associated with a given altitude output measurement is 
limited by the factors given below. 

3.1.1 Timing Uncertainty 

The conversion from AMT to GMT produced a residual 
bias and random error of about 15 ms and -13 ms respectively 
After a correction for the bias error was applied to the SL- 
3 and SL-4 data, there still remained a random error in the 
emphemeris time corresponding to a position error of -100 m. 

3.1.2 Orbit Uncertainty 

Errors in the Skybet orbit are -700 m along-track, 

-270 m cross-track, and -200 m in height [Note 1] . These 
errors produce an uncertainty of about 1 km in position and 
the displacement is mainly along the ground track of the 
spacecraft. 


[lj Skylab EREP PI Data Mtg. at JSC, 16-18 July 1973 


28 


J 


I 


I 


I 


A refined orbit computed by the Geodyne program at Wolf 
Research Corporation for NASA/Wallops is estimated to 
reduce this position error to about -100 m. A comparison 
between the Skybet and Geodyne orbit position values for 
Pass 28 data show a position bias for the Skybet orbit of 
+1500 m with an rms spread of -460. 

3.1.3 Range Tracker Data Smoothing 

Over terrain where the measured heights from indivi- 
dual radar returns may be assumed to vary randomly, the 
smoothing operation of the range tracker on the individual 
range estimates from each radar return produces an addi- 
tional time delay in range readout. For the 2 Hz nominal 
bandwidth of the range tracker, this means an additional 
time delay of 1/4 sec to the time tag for the measured 
height value. This time delay could further increase if the 
antenna angle is off nadir, since the slope of the received 
waveform is reduced, with a corresponding reduction of the 
range tracker bandwidth. Thus, for typical cases, the 
position of the observed measured height as indicated by the 
time tag may be shifted by anywhere between 2 km and 8 km due 
to the range smoothing by the range tracker; Examples of 

this shift will be further discussed in Section 5. 

% 

3.1.4 Radar Footprint 

The position of the reflecting area for which a 
height measurement is obtained may be located anywhere 
within the radar footprint. The uncertainty of this 
position is estimated to be of the order of ±2.0 km, and 
it is present either in -the case of dominantly specular 
returns in flat. areas, or in the case of diffuse returns 
in complex terrain. 
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3.1.5 Antenna Off-Nadir Angle 

The altimeter antenna is normally directed to the 
nadir position. However, in operation the antenna pointing 
is anywhere between 1/2 and 1° from nadir, which can produce 
position shifts between 4 and 8 km. 

A summary of the individual position errors is given 
in Table 7 . 


TABLE 7 


POSITION ERRORS 


Timing 

Orbit 

Range Tracker 
Footprint 
Antenna Angle 


±100 m rms 
±100 m rms 
-4 km ±4 km peak 
±2 km rms 
±8 km oeak 


3.2 Height Error 

The height measured by the range tracker of the radar 
altimeter over terrain is determined from a small portion of 
the radar return where a proper waveform is present. Any 
changes in amplitude, pulse shape, and range will increase 
the range jitter and reduce the precision of the height 
measurement. 

3.2.1 Amplitude Variations 

There are often large rapid amplitude fluctuations 
observed over rough terrain, which cannot be fully controlled 
by the AGC. In the case of sudden drops in signal strength, 
the range tracker may operate for short intervals-on just 
the receiver noise, while it searches for a signal of the 
proper waveform to preserve lock. This search can produce 
peak height errors of the order of 100 m, until the AGC 
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increases receiver gain sufficiently for the detection of 
signal. In some cases where the ramp is displaced, the 
recovered signal will be outside the split-gate of the range 
tracker, and range lock is lost. 

3.2.2 Pulse Shape Variations 

Large changes of pulse shape due to a change in 
reflection mechanism will also increase height jitter. In 
addition, the smaller amplitude and less-defined waveform 
for a diffuse return may further increase range jitter. The 
magnitude of such errors will be similar to those caused by 
amplitude variations. 

3.2.3 Range Rate 

When very rapid and large changes in height cannot be 
followed by the range tracker, large height residuals may 
be produced. These errors, depending on the rate and 
magnitude of the height changes in the topography, may 
produce errors up to 200 m. 

3.3 AGC Error 

The variation in AGC output is proportional to the mean 

of the peak received power. However, the readout values may 

be in error by several JB due to improper calibration of the 

AGC system. An uncertainty of a few dB is inherent in the 

% 

AGC measurement, but since only relative values were used 
in the interpretation of the results, no additional effort 
to reprocess data, including a correction for this error, 
was considered necessary. 

3.4 Errors in Waveform Analysis and Statistical Distribution 

The errors in S&H gate output are mainly due to range- 
tracker jitter, which causes a displacement between the 
return signal and the S&H gate positions and consequently a 
"displaced" sampling of the waveform. The jitter, which 
reflects the. attempt of the range-tracker loop to track the 
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noisy signal will increase as the signal-to-noise ratio 
(SNR) decreases, and the errors in S&H gate output will vary 
accordingly. Individual gate bias levels will also contri- 
bute to gate output error for small SNR. 

' 3.5 Topographic Map Height Errors 

The precision of height measurements from topographic 
maps is usually considered to be about one-half of the 
contour intervals on the map, which varied from —30 to 
-10 m depending on the complexity of the terrain. 




4 . 
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DESCRIPTION OF ALTIMETRY DATA AND TERRAIN 
CHARACTERISTICS 

4.1 General Description of Observations 


Data on a total of 13 Skylab altimeter passes over 
"terrain were received by NRL. A quick survey showed that 6 
of the 13 passes contained useful data, and five of these 
were analyzed in detail. A list of the pass numbers for 
each of the Skylab missions received by NRL is given in 
Table 8. 


SL-2 

SL-3 

SL-3 


TABLE 8 

PI, Pll 

P16 , P17 , P18, P23 , P28, P31 

P4 , P6 , Pll, P37, P49 


. The 6 passes which acquired relatively reliable terrain 
data are underlined. Pass 11 of SL-2 and Pass 31 of SL-3 
contained only calibration data. All SL-4 passes were 
marginal due to a large drop in received power (-24 dB) 
caused by apparent antenna gain deterioration. 

The most reliable operation was achieved in passes 16, 

17 and 18 of SL-3 where all 3 sub-modes were activated. For 
the remaining passes (1 of SL-2, and 23 and 28 of SL-3) only 
the first 2 sub-modes were employed, due to a large number of 
frames where range tracker loss occurred. Table 9 provides 
a summary of the number of frames acquired in each pass and 
the approximate location of the ground track for each sub- 
mode. The ground tracks of the 5 passes which were processed 
by NRL are shown in Figure 12 and a general description of 
the terrain observed on each pass is -given in -the. following 
paragraphs. 
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Pass 1 (SL-2 ) 5/30/73 


The first pass of SL-2, travelling southeastward, 
acquired range lock over the Pacific Ocean and crossed the 
West Coast of the United States at Coos Bay, Oregon. The 
ground track traversed a variety of terrain with complex 
topography, which included large mountain ranges, forests, 
canyons, and some valleys and small mountain lakes. The 
Observations covered portions of Oregon, California, Nevada, 
and Utah and were terminated over the Marble Canyons in 
Arizona. As indicated in Table 9, only 25 percent of the 
available frames (data blocks) had range lock. This pass is 
typical of the Skylab altimeter performance over complex 
terrain, where a narrow pulse range tracker could not 
accommodate large changes in both topographic height and 
received power. 

Pass 16 (SL-3) 9/6/74 

The observations of Pass 16 started over Baja 
California in Mexico, and continued over the Gulf of 
California, and the Desierto del Altar. The ground track 
then entered the United States near Ajo in Arizona, where it 
covered some valleys and mountains. It then continued over 
the highlands of Arizona and New Mexico which are mostly 
covered by ranchlands, and over large mountains in Colorado. 
The pass terminated over a flat area near Pueblo, Colorado. 
This was the best run as a result of a combination of reliable 
instrument performance and gradually changing terrain 
features. More than 90 percent of the frames had range 
lock. 

Pass 17 (SL-3) 9/7/73 

The ground track of Pass 17 traversed mostly smooth 
areas and in spite of poor antenna pointing range lock was 
achieved for more than 90 percent of the available frames. 
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The outstanding terrain features included farmlands in Iowa, 
swamps in Wisconsin, and two large lakes. Lake Michigan and 
Lake Huron. 

Pass 18 (SL-3) 9/9/73 

The ground track of Pass 18 is similar to that of Pass 
17 in terms of relief. Large areas of flats and ranchlands 
in Texas, gently sloping hills in Oklahoma, Arkansas and 
Missouri are typical terrain features. There was range lock 
on about 90 percent of the frames. 

Pass 23 (SL-3) 9/11/73 

The Skylab altimeter in this pass was operated over 
Idaho, Montana and North Dakota. About 75 percent of the 
pass had range lock. 

Pass 28 (SL-3) 9/13/73 

The ground track of Pass 28 includes a combination of 
the type of terrain observed in Pass 1 of SL-2 and Passes 
16, 17 and 18. The run started in the mountains and forests 
of Nevada, then passed over the Salt Flats and desert areas 
of Utah, and finished over a mixture of mountains, plateaus 
and ridges in Idaho and Wyoming. About 50 percent of 
available frames had range lock. 

4.2 Fine-Structure AGC and Ground height 

The fine-structure Skylab altimeter AGC output, as 
defined in paragraph 2.2, represents the combined effects of 
complexity of terrain and reflectivity along the ground 
track. Low relief, valley floors, and specular reflection 
will enhance AGC signal level, while high complex relief and 
low reflectivity or signal scattering terrain ‘types,’ inclu- 
ding vegetation, will produce low AGC signal levels. In 
other words, AGC signal level will generally vary with the 
particular type of terrain in this prescribed manner. Low 
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signal levels will also increase range noise, and in extreme 
cases will result in loss-of-lock by the range tracker. 
Typical examples of such relations for various types of 
terrain are shown in Figures 13 through 19. 

Pass 16 

In Figure 13 is shown the altimeter output for the 
first part of Pass 16. The upper curve is the AGC response 
to the observed terrain, and the lower curve is the corre- 
sponding altimeter tracking profile. The ground track, 
starting in the Valle Chico in Mexico (DB 180) , passes over 
the Sierra San Felipe mountains (DB 183) , and then descends 
to the Gulf of California (DB 184, 185). The spacecraft is 
over the Gulf of California for about nine seconds (DB 186- 
195) and it then remains over the Desierto del Altar for the 
remainder of the plot (to DB 204). It is apparent that the 
presence of the San Felipe mountains reduces the receiver 
power by about 25 dB. The initial sudden drop in AGC signal 
level results in a loss of range lock for one frame (DB 
181), but because of good antenna pointing, range is 
recovered in DB 182, and the cracker even follows the rela- 
tively large descent of about 500 meters to the coastal 
flats of the Gulf. The same altimeter behavior can. also be 
better observed in Figure 22 which includes a topographic 
map profile for the same ground track. 

Over the Gulf, the AGC signal level increases to and 
remains a constant value of -30 dB until the Desierto del 
Altar is reached. The altimeter-measured ground height 
remains stable over the Gulf and then increases smoothly 
with the topographic map land elevations as can be seen in 
Figure 22. 

The above case is a good example of the potential 
performance of the instrument over land and water areas of 
moderate relief. 
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Pass 17 

Another example of fine- structure altimeter output is 
shown in Figure 14 for Pass 17 over Wisconsin. The terrain 
consists mostly of swamps and small lakes, but here there is 
a large fluctuation of the AGC signal of 10 to 20 dB, which, 
together with the lower signal level, is due to a large off- 
nadir antenna pointing angle as will be explained later m 
Section 4.4. A correspondence between the AGC signal level 
and the altimetry profile is less evident in this case. 
However, strong signals usually reflect flat or valley 
terrain, while a drop in signal power means steeper relief, 
with the magnitude of the fluctuation being enhanced by poor 
antenna pointing. In spite of poor performance at such 
angles, the range tracker remained locked in over the 
relatively even terrain, and range jitter, while larger than 
that in Pass 16, was still relatively small (Note increase 
of elevation scale by a factor of 4.) 

Another sample of Pass 17 over the Lake Michigan region 
is shown in Figure 15. In this case, the spacecraft 
approached the lake from the west over farmland in Wisconsin. 
The altimeter profiles the approach to the lake, but the 
high frequency variations in height represent range tracker 
jitter due to high range rates. The high AGC values over 
the western portion of the lake and the later decrease are 
probably due to an initially smooth water surface which is 
subsequently disturbed by prevailing westerly winds as the 
eastern coast of the lake is approached (NOAA Daily Weather 
Map-i, 3-9 September 1973) . 

In Figure 16 the altimeter output for the Lake Huron 
region is shown. Here, again the altimeter follows the 
descending elevation to the lake. However, in this case, 
the AGC signal level over the lake is about 15 dB below that 
of Lake Michigan. This is produced by a sub-mode change 
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which tilts the antenna by about 1/2 degree in pitch, but as 
will be shown later, the pointing is actually off by 1 
degree due to an initial offset in the pointing of the 
antenna beam. The large off-nadir angle also results in 
large range jitter over the lake region. 

Pass 28 

Other examples of fine-structure altimeter outputs are 
shown in Figures 17 through 20, where cloudless Skylab 
photography taken over the ground track is also included. 

In Figure 17 the terrain consists of alternating mountain 
ranges and valleys in Nevada and several range— lock losses 
are observed. The generally low AGC signal and the large 
change in AGC over the Cherry Creek Mountains are apparent. 
It is seen that even for a constant AGC signal level, the 
range acquisition is delayed if the value falls below -50 
dBm. The run is continued in Figure 18 and reliable opera- 
tion is attained only after large AGC signal is regained 
over the Bonneville Salt Flats in Utah. 

Figure 19 again shows the potential performance of the 
altimeter for large and constant AGC. The first drop in AGC 
due to the Newfoundland mountains has no adverse effect on 
range measurement since the change in AGC is less than 5 dB 
and the AGC signal is quite strong. A later drop in AGC 
signal over the lakeside mountains near the Great Salt Lake 
is much larger (20 dB) and there is a corresponding large 
jitter in measured height before the AGC and the range 
tracker restabilized. In Figure 20 an increase in jitter as 
well as occasional range tracker loss are indicated as the 
Skylab passes over complex relief in Wyoming. Reliable 
operation is again achieved over flats and the Fork River. 
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4.3 


Comparison of Radar Altimetry with Map Topography 


Altimeter-measured and topographic map elevations along 
the. ground .track- were- compared by- the -following- -methodr 


which was finally adopted after a few trial and error 
attempts. This consisted of estimating maximum and minimum 
map elevations in a footprint in order to establish bounds, 
and then comparing altimeter-determined ground heights with 
these bounds. Height residual is defined here as the value 
by which the altimetry ground height exceeds these bounds. 

The height residual could then be positive or negative. The 
relation of the altimeter-determined elevations to the 
maximum-minimum bounds is depicted in Figure 21 through 25, 
while the height residuals are shown in Figures 26 through 
30 and also tabulated in Tables A-l through A-5 for each data 
block. In this instance, fine- structure altimetry data was 
averaged over a data block or 1.04 seconds in time to reduce 
AGC and range jitter, and then compared with the map data 
for the corresponding block. 


This method of analysis showed that large height 
residuals • existed for several data blocks prior to and 
following range lock loss and also when large off-nadir 
antenna angles existed. In retrospect, this technique 

S 

provided a general means of evaluating altimeter performance 
in terms of accuracy, response and reliability. The follow- 
ing sections show the results of this operation and the 
possible causes of any anomalies observed. 

Pass 1 

In Figure 21 the topographic heights for Pass 1 are 
shown plotted vs. time. As has been indicated previously, 
this pass scanned a region of very complex relief and range 
lock was lost over many areas. In general ,' however , the 
altimeter-measured heights closely followed the topographic 
map profile. There are a few major discrepancies which will 
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be discussed later. It is seen that range lock is lost 
initially as the Skylab ground track crosses the West Coast 
of the United States at Coos Bay, Oregon. However, as the 
relief becomes more uniform, range is again acquired, but it 
is lost again as the height increases rapidly . Only for 
more or less uniform relief can the range tracker follow the 
height changes. Large height residuals are often observed 
near large changes in height, and the range tracker tends to 
remain inertially locked— in to a previous elevation. This 
can be clearly seen near Alturas, California where altimeter- 
measured height is radically below map height for _he case 
of rapidly rising elevation. Here, the altimeter-measured 
heights correspond to the lower map elevations within a 
footprint. Contrarily, when the height is decreasing as 
near Cedar City, Utah, the range tracker tends to lock into 
the preceding, in this case, higher elevation. This action 
is equivalent to a low pass filter in the system which 
produces a slight lag of about one footprint in altimeter 
height measurement for rapid range rates. 

Pass 16 (Figure 22) 

The performance of the altimeter in Pass 16 shows the 
high accuracy and reliability which can be attained for 
terrain of moderate relief in conjunction with accurate 
nadir pointing. The decrease in height as the Gulf of 
California is approached as well as the gradual increase in 
height of the Desierto del Altar area are clearly seen. The 
tracking of the plateau region in Arizona is evident and 
range lock is lost only when a rapid and large increase in 
height at Mesa, Arizona occurs. Reliable operation is 
obtained over the highlands in New Mexico, but the steep 
relief near Durango, Colorado cannot be followed by :he 
altimeter and range is lost for a considerable distance. In 
spite of the occasional range loss, this pass represents the 
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best of the height observations made by the Sky lab and it 
also had the smallest rms error. 

Pass 17 (Figure 23) 

The ground track of this pass covered low areas of low 
relief and operation was relatively reliable in spite of a 
known off-nadir pointing angle. The early part of the pass 
was over farmlands in Iowa and Wisconsin and there is good 
correlation with map topography. Some residual errors are 
observed, however, near large vertical structures such as 
those near Dubuque, Iowa and Lacrosse, Wisconsin. The lag 
effect in altimeter ranging is again evident near the 
decrease in height as Lake Winnebago, Lake Michigan and Lake 
Huron are approached. The constant bias over Lake Huron is 
due to off-nadir pointing and will be more fully explained 
in the next section on height residuals. 

Pass 18 (Figure 24) 

Pass 18 is similar to Pass 17 in terms of observed 
topography and again the altimeter height profile follows 
the map elevations. Surprisingly, there was reliable 
tracking of the large height changes over Texas, but this is 
because of moderate relief and range rate manageable by the 
altimeter. Range lock is lost near ^Russellville, Arkansas, 
where steep relief persists over several frames. Small 
deviations near Ardmore, Oklahoma are not due to height 
changes but are caused by a sudden decrease in reflectivity 
as will be shown later in Figure 29. 

Pass 28 (Figure 24) 

The ground track of Pass 28 covers complex terrain in 
the early and late portions of the observations and homo- 
geneous terrain in the middle portion. Range lock is lost 
over mountains near Ogden, Utah, and recovered reliably only 
in a dry basin near Landon, Wyoming. Major discrepancies 
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were observed near Ely, Nevada, due to a "lag" effect, and 
near Great Salt Lake and Preston, Idaho where large AGC 
changes were recorded. In general, the altimeter again 
measures the lowest elevations within the observed footprint. 

4.4 AGC and Height Residuals 

In Figure 26 through 30 are shown smoothed AGC values 
and associated minimum height residuals. The upper curve in 
Figure 26 represents the AGC values for Pass 1. It is seen 
that, in general, AGC values are below -35 dBm except over 
Klamath Falls, Oregon, and the Black Rock Desert, Nevada. 

The lowest values are from forests and mountains, with 
intermediate values for valleys and hills. The AGC values 
also .show considerable fluctuations even after smoothing. 

This produces frequent range lock loss as seen in the lower 
curve. The larger negative residuals are always associated 
with a large drop in AGC values as seen for the plot near 
the Madoc National Forest, also after the Black Rock Desert, 
and for the valleys and hills in Nevada (DB 220) . The AGC 
value for the Pacific Ocean (-35 dBm) , shown at the beginning 
of the plot, is slightly lower than the average for the 
ocean indicating some off-nadir pointing, since relatively 
smooth sea-states were present. 

In Figure 27 are given the results for Pass 16. The 
fluctuations of the AGC values are considerably reduced and 
a reference value of -30 dBm for the Gulf of California is 
obtained. While mountains and forests still produce small 
AGC values (-50 dBm) , the small fluctuations in level allow 
the range tracker to follow the changes in height without 
range lock loss. The few failures in range tracker opera- 
tion include one due to a drop in AGC (DB 242) and one due 
to a large vertical structure in DB 315 as indicated pre- 
viously. The large height residuals over ranchland (DB 
270 - 290) are caused by a combination of low AGC signal 
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levels and 100 MHz bandwidth operation in sub-mode 1. 

The height residuals and the corresponding AGC values 

_f or pass -17 are— shown-in— Figure— 28,— -The-residual- values - — 

have been generally greatly reduced, but again some of the 
larger negative residuals are associated with a decrease in 
received power. The case of the large positive residual 
(45 m) is probably due to a combination of the large elev- 
ation in the adjacent footprint and the low AGC value, which 
permits the weaker but earlier signal received from an off- 
nadir angle to dominate the range tracker operation. The 
low AGC causes the residual to reverse direction in the next 
frame which is then followed by a loss in range tracker 
lock. Another significant height difference is found over 
Lake Huron where a height bias of about -10 m relative to 
Lake Michigan is indicated. According to the map, the two 
lakes should be at the same elevation. investigating the 
performance of the altimeter during this interval, it was 
discovered that (1) the AGC level for Lake Huron is about 
-15 dB below that for Lake Michigan, and (2), a sub-mode 
change occurred over Michigan that shifted the antenna pitch 
angle by 1/2 degree. If it is assumed that before change of 
sub-mode, the antenna angle was offset from nadir, a situa- 
tion that occurred frequently over ' the ocean and which was 
later confirmed from waveform analysis over Lake Michigan, 
and also that the sub-mode change increased the off-nadir 
antenna angle to one degree, then an increase of about 10 m 
in satellite height would occur. This bias is obtained from 
a plot of height bias vs. antenna angle in a G.E. Report 
Skylab Altimeter , March 19, 1971, Figure 7, p. III-52. The 
initial off-nadir angle also accounts for the large AGC - 
jitter over the farm areas in Iowa. 

The AGC values and height residuals for Pass 18 are 
shown in Figure 29. The AGC jitter is greatly reduced due 
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to better nadir pointing and less relief. High AGC values 
indicate water surfaces over Lake Texoma, the White River in 
Missouri and the Mississippi River in Illinois. The Ozark 
Mountains in Missouri, forests in Arkansas and small mountains 
in the 'Northern part of Texas coincide with low AGC values, 
and the ranchlands in Texas provide intermediate values of 
about -35 dBm. The height residuals are primarily due to 
rapid height changes which cause preceding height values to 
be measured by the range tracker. 

The values of the height residuals compared with the 
AGC values for Pass 28 in Figure 30 again show the correla- 
tion between large negative height residuals and received 
power. The large negative height residuals are almost 
always followed by a loss of range track. Figure 30 shows 
the considerable AGC fluctuations that occur over complex 
terrain. The large variations are partly caused by signifi- 
cant contributions from specular reflections which produce 
strong signal returns over short intervals and which just as 
quickly disappear. Typical of the large returns that are 
possible from specular reflections is the received power for 
the Great Salt Desert which saturated the receiver. 

4.5 Waveform Analysis 

The results of waveform analysis are described in 
Figures 31 through 41. The first two figures give the 
expected and measured waveform for returns from the ocean 
surface as a function of off-nadir angle. The other diagrams 
give plots of signal waveforms for different types of 
terrain, averaged for 100 pulses, and also the normalized 
auto-correlation functions for 100 samples of No. 4 
sample-and-hold gate output. The jump or lag index in the 
latter is in units of 0.01 seconds of time, which is 
equivalent to 72 meters in ground track distance. The 
criterion for correlation distance or the size of a land 
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patch which gives a coherent signal return is a normalized 
correlation value larger than about 0.2 to 0.25. A steep 

correlation curve dropping to zero indicates a short corre- 

lation distance and a randomly scattered signal. On the 
other hand, a stretched out curve indicates a long correlation 
distance and large specular reflecting areas. 

The theoretical ocean radar return consists of a 
triangular pulse whose rise time is about 100 ns and which 
decreases gradually to about 20 percent of the maximum value 
over the next 30C ns. As the off-nadir angle increases, 
there is a small decrease of rise time, but the trailing 
edge decreases to only about 35 percent of maximum for a 0.5 
degree and 75 percent of maximum for a 0.75 degree off -nadir 
angle. The measured waveform for the three sub-modes of 
mode 1 for the ocean surface (Pass 11) shows that (1) there 
appears to be an off— nadir angle even for sub-modes 0 and 1 
which is about 0.5 degree, (2) as sub-mode 2 is initiated, 
the 0.5 degree pitch angle further increases off-nadir 
operation to about 0.75 degrees, as determined from the 
trailing edge, (3) there is a shift in the range tracker 
position as is apparent from the shift in the rise time for 
0.5 degree pitch angle change, and (4) there are considerable 
fluctuations in the radar return particularly in the trailing 
edge of the waveform, which is due to intrinsic signal noise 
since only 100 radar returns were averaged. These waveforms 
are used as a reference for comparison with waveforms 
obtained from terrain reflections. 

Figure 33 shows the noise- and ramp-position gate 
outputs of Pass 1 for two adjacent frames, as sub-sub-mode 0 
is changed to 1, over a relatively complex mountainous 
terrain in Nevada. When compared to an ocean return, it is 
apparent that there is some structure in the rise time which 
consists of both scatter and large-body reflection. The 


large-body reflection is indicated by the spatial auto- 
correlation function shown on the right which shows that 
horizontal smooth areas of the order of about 150 m are 
present. 

Similar radar returns are shown in Figures 33 and 34 
for Pass 16. Note that noise- and ramp-gate outputs are 
shown as before in Figure 34, but that ramp- and plateau- 
gate outputs are shown in Figure 35. 

The radar return in Figure 34 is from the Gulf of 
California and the rapid decrease in signal strength after 
the peak indicates reliable nadir pointing, which was also 
observed from other data for this pass. Furthermore, the 
auto-correlation curve shows that pure scattering dominates . 

4 . ... . 

Figure 35 shows two returns for the 10 and 100 MHz 
submodes, which respectively reflect from mountains and 
ranchland in Arizona. In this case, the valleys in the 
mountains give a dominantly specular return as indicated 
by the rapid drop-off of the plateau region of the return 
waveform and the correlation distance of about 100 m. On 
the other hand, the ranchland appears to be a more diffuse 
target. 

In Figure 36 the high correlation indicates that areas 
of over several hundred meters length contribute to specular 
return in passing over farmland in Iowa and swamps in 
Wisconsin. Another waveform for returns from farmland in 
Iowa is compared with the dominantly diffuse return from 
Lake Michigan in Figure 37. The high trailing edge of the 
Lake Michigan return indicates a large off-nadir angle as 
was indicated previously for Pass 17 observations. Another 
difference between a terrain and a water return is the much 
larger fluctuations in the trailing edge of the water return 
compared to large-body reflection from terrain. 


Waveforms for Pass 18 observations are shown in Figures 
38 and 39. While in Figure 38 both returns are from 
ranchlands in Texas, the presence of brushwood in one area 
has considerably increased the scattering component. Large 
body reflections of at least 300 m extent are indicated in 
the other return. In Figure 39 a strong specular return is 
noted over flats in Oklahoma. 

The wavefrom results of Pass 28 are given in Figures 
40 and 41. The mountain return in Figure 40 shows both 
scattering and large body reflections, with scattering 
apparently the dominant mechanism. In the following figure, 
the pulse shape indicates an almost purely specular return 
from the salt flats in Utah but the relatively low auto- 
correlation value also seems to show that large area 
reflections can fluctuate rapidly. 

4.6 Catalogue of Terrain Signatures 

The radar signatures of the Skylab ground tracks for 
Passes 1, 16, 17, 18, and 28 are listed in Tables A-l through 
A-5 in the Appendix. Here, the AGC signal level, the 
reduced altimetry ground height, and some waveform analysis 
data are tabulated against topographic map elevations 
accompanied by a brief geographical description. An expla- 
nation of the column headings is given in the Appendix. 

4.7 Available Data Base 

For the reference of those who may be interested in 
further analysis of the data or a more detailed description 
of a particular area, a set of five volumes is available at 
NRL that contains all the computations performed by NRL on 
the five Skylab altimeter passes. Each volume is divided 
into three parts. The first part contains the initial 
evaluation data described previously in Section 2. The 
second part lists the following outputs for each data block: 
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(a) 4 AGC values and mean 

(b) 8 measured range values 

(c) 8 corrected range values and mean 

(d) 8 corrected spacecraft heights and mean relative 
to geoid 

(e) 8 height residuals and mean 

(f) Geoid height. 

In the third part, the S&H gate outputs and related 
statistical analyses are recorded. For each data block and 
S&H gate output, the following information is available: 

(a) 104 S&H gate amplitudes after correction 

(b) the mean and standard deviation of these amplitudes 

(c) the amplitude distribution 

(d) the auto-correlation of the amplitudes. 
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5. DISCUSSION OF RESULTS AND CONCLUSIONS 

5.1. Radar Return Characteristics of Terrain at Normal 
Incidence 

An evaluation of the Skylab radar altimeter terrain 
data indicates that certain radar return characteristics are 
associated with certain geometrical and physical properties 
of the observed terrain. The more significant of these 
characteristics and relations for normal incidence obtained 
from the Skylab altimeter observations are discussed below. 

5.1.1 Equivalent Radar Cross-section of Terrain at Normal 
Incidence 

For the radar altimetry data used in this particular 
analysis, a distribution of radar cross-sections such as 
shown in Figure 42 was obtained. This distribution simply 
reflects the relative occurrence of the different types of 
terrain sampled and it is therefore a signature of regional 
geography or terrain types along a particular ground track 
or tracks. On the other hand, however, this distribution 
and the data in Table 10 are indicative of the range of 
cross-sections to be expected over mixed land-and-water 
areas. It should be noted that the loss of range-tracker 
lock and data over steep mountainous relief with vegetative 
cover which often have small cross-sections of the order 
of 60 dB- above-one- square-meter (dBsm) , account partly for 
the small number of data samples taken over such terrain, 
whereas in reality, mountainous terrain could be the domi- 
nant terrain type in some region. This means that careful 
interpretation of radar cross-section signatures is neces- 
sary. 

The actual equivalent normalized radar cross-section 
for a terrain type may be considerably larger than shown 
since the effective scattering area at a given height is 
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smaller than the footprint in areas of complex relief. The 
upper and lower bounds of the radar cross-section are 
limited by the saturation and noise levels of the receiver 
respectively, but an extrapolation of the data base indicates 
that about 70 percent of the terrain radar returns would be 
included in the distribution. The saturation level of the 
receiver, that is, 110 dBsm, was only exceeded over the salt 
flats in Utah. The relatively large radar cross section 
over many complex terrain areas confirms the results of the 
waveform analysis, which had indicated that the radar return 
was produced by a dominantly smooth-area reflection at 
normal incidence, even at the short wavelength of 2.15 cm. 
While the indicated values of the radar cross section 
correspond to the maximum returns for the observed areas due 
to the power-level selecting mechanism of the AGC loop in 
the radar receiver, there may be additional significant 
power contributions from other height levels within an 
observed footprint, which signals are outside of the AGC 
gate and therefore not recorded. 

5.1.2 Relation of Received Power to Terrain Features 

The relative received power as a function of the 
observed geometric and physical terrain features is listed 
in Table 10. It is seen that only in' flat desert areas and 
for smooth water surface does the received power exceed that 
from ocean return, while for ridges and dry lakes, it is 
about 10 dB below that for the ocean. The presence of 
vegetation, such as farmland or brushwood reduces this value 
by about 5 dB. The relatively lower values for areas with 
vegetation indicates that the vegetation reduces the specular 
component of the radar return. The smallest values of 
received power are obtained over mountains and forests. 

There is considerable scatter in these values since the 
received power varies with both terrain geometry, soil 
moisture, and vegetation. 
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TABLE 10 


Dynamic Range of Peak Terrain Radar Return 
Relative to Ocean Return 


Salt Flats 
Lakes 

Ocean, Deserts 

Valleys, Plains, Cities, Swamps 
Ridges, Canyons, Dry Lakes 
Hills, Ranges, Mountains 
Cliffs, Forest 


dB 

20 

10 

0 

- 5 
-10 
-12 
-16 


Although off-nadir antenna pointing degrades both the 
power and the quality of the return signal, the effect of 
the pointing angle on return power can be generally quantrfie 
only in the case of returns from continuously uniform surfaces 
such as posed by bodies of water (Section 4.5 and Figures 31 
and 32). Over terrain, the effect is less clearly defined 
because of the combined influence of relief and reflectivity. 

5.1.3 Range Measurement Characteristics 

It has already been indicated in Section 4 that the 
height measured by the Skylab altimeter corresponds to the 
bright spots of the illuminated footprint along the space- 
craft ground track. These bright spots are usually found at 
heights which have the highest probability of normal reflec 
ting areas. The bright spots can range in size from those 
of about 100 meters, corresponding to valley bottoms in 
mountainous terrain (Figure 33) to larger areas of a few 
hundred meters dimension in grasslands, swamplands , and 
similar uniform flat expanses (Figure 36). From this 
standpoint, the radar profiles the height variations of 
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generally the lowest areas within a footprint, and the finer 
undulations of the terrain are smoothed out. The altimeter 
therefore acts as a low-pass filter whose bandwidth is 
considerably smaller than that obtained from the basic 
footprint. 

To show this more clearly, a section of Pass 16 has 
been analyzed in more detail. Figure 43 shows the profile 
obtained from topographic maps along the spacecraft ground 
track with elevations read every one-fourth of a footprint 
step. The bottom curve in Figure 44 represents a five-point 
running average of the former, which filters out wavelengths 
shorter than one-half footprint. When this is compared with 
the middle curve, which represents the height measured by 
the altimeter, it is seen that the high frequency content of 
the altimetry profile has been effectively suppressed in the 
radar measurements. Correlation studies of the two profiles 
indicate that the low frequency components of both profiles 
produce a relatively large correlation of 0.63. However, 
when the spacecraft position is shifted by about -7 km along 
the ground track, the correlation increases to 0.92. That 
such a shift may be necessary is further born out by analyzing 
the differential heights between the mean topographic profile 
and the altimeter measurements. Foi the same shift, the 
mean difference is decreased from 33 m to 17 m and the 
spread in the height difference is reduced from ± 50 to 
±38 m. These results are plotted in Figure 45. This 
apparent shift would indicate that either a time or position 
error in the ephemeris, which was not corrected, exists. A 
part of the shift is due to the range-tracker time-delay lag 
discussed in Section 2. To check the consistency of this 
result, a second section of P16 was analyzed similarly, and 
the results are shown in Figures 46 and 47. In this case 
the average profile from the topographic maps was relatively 
smooth, and the altimetry profile almost coincided with the 
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former. The correlation in this case increased from 0.94 to 
0.98 and the mean difference and spread decreased from 27 m 
to 2 m, and +51 m to ±7 m respectively for the same shift. 

Thus it appears that if the topographic map height is known, 
the position of the satellite over complex terrain could be 
obtain-'-'! independently with an accuracy of the order of 100 
m, by comparing short sections of the map profile over the 
ground track. An analysis of the gross characteristics of 
the differential height given in Section 4 showed that a 
shift in position of -7 km reduced the height differences 
and spread considerably also for P16, 17, and 18. The 
results are summarized in Table 11. 

That the high frequency components of the altimetry 
profile affect the power measurements can be seen by comparing 
the top and bottom curves of Figures 44 and 46. However as 
previously mentioned in Section 2.3.5, no simple relation 
between terrain geometry and power variation could be estab- 
lished, since the latter is affected by a combination of 
factors such as relief, soil moisture, vegetation, and 
surface roughness. 

5.1.4 Waveform and Statistical Characteristics of Radar 

It has been shown in Section 4 that the pulse shape of 
the radar return indicates a dominant large body reflection. 
The drop-off of the correlation as shown in Figure 38 
indicates that the flat areas are of the order of a few 
hundred meters length. This is shown more clearly in Figure 
48 where the amplitude for each pulse return is plotted for 
different sample-and-hold gates in the ramp region. The 
coherent return over several pulses is evident, and the 
fluctuations in amplitude are due to the varying size of 
smooth normal areas. The second half of the plot shows the 
corresponding plateau region -vnd the absence of a signal 
implies that no scattering effects are observed. 
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Previous terrain radar observations from aircraft 
platforms near nadir indicate that the radar return even at 
lower frequencies (480 MHz, 3800 MHz) appears to be predomi- 
nantly diffuse (Reference 9). The dominant specular return 
from orbital heights may be due to the large platform height 

and small pulse width which emphasize the coherent component. 
It has been shown (Reference 10) that if smooth patches of the 

order of 100 m exist, the ratio between the specular and 

diffuse component is given by 



where , k = the reflection coefficient, 

h = the height of the observation platform , 

a = normalized radar cross section , 
o 

c = velocity of light, 
and T = pulse width. 

Thus as the platform height increases, and the radar pulse width 
decreases, the specular component will become dominant. 

When the coherent return is reduced sufficiently, the 
range tracker loses lock. An example is given in Figure 49 
for DB 243 of H6 where the amplitude in ramp gate 7 is 
shown. The drop in amplitude is due to the loss of signifi- 
cant signal from smooth areas due to rapidly increasing 
height and large increases in slope. Before the range 
tracker loses lock, the altimetry ground height is less than 
the map height, since the tracker persists in searching near 
the same height level. 

5.2 Radar Altimeter Performance 

While the radar altimeter lost range lock about 30 
percent of the time in the five analyzed passes, the remain- 
ing 70 percent of the data showed reliable altimeter 
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performance over many complex terrain features. The flexibi- 
lity and high sensitivity of the Skylab altimeter made it 
possible to obtain height data that reasonably closely 
profiled the low frequency characteristics of the relief 
along the spacecraft ground track as already discussed in 
Section 5.1.3. The limitations to its performance were 
primarily due to the range tracker and the AGC loop responses 
and interaction between the two. A combination of off-nadir 
antenna pointing and the time-delay truncation of the radar 
return by the range tracker and AGC processor degraded 
return signal amplitude and quality. The situation worsened 
with increasing relief of the observed terrain, and in 
extreme cases, resulted in loss of range locks. 

The off-nadir pointing angle of the antenna can be 
obtained from the radar return characteristics from water 
surfaces by measuring either the rate of drop-off of the 
trailing edge in the plateau region or the amplitude and 
bandwidth of the range jitter. Since normally over terrain 
only small water surfaces are present, the signal plateau 
region usually could not be used. An exception was the 
observation over Lake Michigan in Pass 17, and as has 
already been indicated in Section 4, the drop-off of the 
trailing edge indicated that the antenna pointed off-nadir 
by about 0.75 degrees. The alternate method (Ref. 11) of 
height- jitter characteristics can be used, if the plateau 
region for signal returns from large water surfaces is not 
available. While this method is only strictly valid for 
large off-nadir angles (> 1°) , it has been applied to a 
sample of each pass for estimating system performance. The 
results are gi\en in Table 12. It is seen that the largest 
bandwidth is obtained for the Gulf of California, which 
implies an off-nadir angle of less than 0.5 degree. The 
reduced bandwidth over Lake Huron of 1.8 Hz would indicate 
an off-nadir angle of about one degree. This agrees with 


59 


1 




other factors, such as the low received power from Lake 
Huron as compared with that from Lake Michigan and the Gulf 
of California, and it is further confirmation of the height 
error of about 10 m discussed previously in Section 4.4. 
Similar analysis, applied to homogeneous terrain areas, are 
shown for small sections of ranchland in Texas (Pass 18) and 
the Salt Flats in Utah (Pass 28). The improvement in jitter 
(0.8 m) for large specular return is apparent over the Salt 
Flats. 

In spite of the above limitations, the Sky lab altimeter 
accommodated large range and AGC rates and became inoperative 
only over the larger vertical height dispersions. 

TABLE 12 


Range Tracker Jitter 

and Bandwidth 



Jitter 

Bandwidth 


(m) 

(Hz) 

PI, Pacific Ocean 

2.9 

2.0 

P 16, Gulf of California 

s 2 

2.3 

P 17, Lake Michigan 

2.8 

2.3 

Lake Huron 

3.25 

1.8 

P 18 Ranchland, Texas 

1.9 

2.2 

P 28, Salt Flats, Utah 

.8 

2.3 


5.3 Design Considerations for a Dedicated Terrain Altimeter 

Terrain observations at nadir require a radar altimeter 
system different from one used exclusively for ocean observa- 
tions. Significant differences in the nature of the target 
are: 
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1. The terrain represents a complex radar target where 
range variations within a practical-size footprint may- 
amount to several hundred meters . 

2. Range fluctuations are rapid and frequent within a 
footprint, and their pattern can vary from one footprint to 
another. 

3. Reflection will be both specular and diffuse, with 
the former being favored since it returns a stronger coherent 
signal. 

4. A number of factors, such as slope, rate of slope 
change (sharp or gentle relief) , dielectric constant, 
surface roughness, and vegetation, afiect terrain reflecti- 
vity, which is in contrast to a relatively simple scattering 
mechanism for ocean surface returns. 

5. Range and reflected power variations within a 
footprint, due to the causes described, could be of the 
order of 1,000 m and 30 dB respectively. 

A full description of the radar return would consist of 
a three-dimensional representation such as shown in Figure 
50, where received power could be shown as a function of 
height z, ground track distance x, and track displacement 
y. The signal sampled by a specific narrow range gate would 
yield the reflected power for a particular terrain elevation. 
This would provide information on both the relief and the 
radar reflection characteristics for each elevation of the 
terrain along the sub-satellite ground track. By combining 
a series of measurements along parallel tracks, a two- 
dimensional reflectivity profile for each elevation is 
obtained. 

The design parameters required of a satellite altimeter 
which would provide such an output are: 

1. A maximum sampling rate which will yield independent 
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data samples. This would make it possible to identify and 
to describe the extent of specular-return regions and to 
locate reflecting areas with height and horizontal resolu- 
tions of the order of one and 100 meters respectively for 
satellite altitudes. 

2. The use of bearawidth- limited operation, by eix-her 
decreasing antenna size or increasing pulsewidth, so that a 
footprint can be uniquely associated with a given land area. 

3. A threshold range tracker which detects the earliest 
return and responds to range rates of about 1,000 m/s. 

4. An adaptive S&H gate with response time of the 
order of one second, where the width of the operating 
sampling gate can be varied to accommodate the maximum 
height range of the observed terrain. For this purpose, 
about 20 S&H gates would be needed to sample the radar 

5. A dynamic IF range of 30 dB, so that both small 
signals from small scattering areas and large signals from 
specular or larger scattering areas can be measured 
simultaneously . 

6. Refined antenna-pointing calibration, which can be 
obtained from the radar return over the ocean or other water 
surfaces for which reflection or scattering characteristics 

can be assumed to be known. 

A possible set of system parameters required to produce 
such a performance for a minimum radar cross section of 
dBsm (200~ns pulsewidth) is given in Table 13. 

A 'date-processing scheme for such a system is shown in 
Figure 51. Three parallel detector circuits, different in _ 
power levels through lOdB attenuator inputs, cover an 
instantaneous dynamic signal power range of 30 dB. The AGC 
is controlled by the peak signal as in the Skylab altimeter, 
but the first detector is also allowed to sense signals 30- 
dB below the peak level. 
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TABLE 13 

System parameters for a dedicated terrain satellite radar 
altimeter with a satellite height h g = 4 x 10 5 m. 
Transmitter 


TWT power output 

2kW 

Wavelength 

2 cm 

Pulse-width 

200 ns 

Pulse Repetition Frequency 
with burst of 10 pulses 

separated by 200 ys 

100 

Average transmitted power 

0.4 W 

System loss 

6 dB 

Antenna 

Paraboloid diameter 

1.5m 

Aperture efficiency 

50% 

Gain 

43 dB 

3 dB beamwidth 

1° 

Receiver 

System temperature 

500°K 

Bandwidth 

5 MHz 

SNR (a = 53 dBsm) 

10 

Data Processing 

Type of acquisition and track 

Threshold 

Tracking bandwidth 

10 Hz 

t (Tracking used to position S&H 

gate #5) 

Number of S&H gates 

20 


S&H gate #1 precedes S&H gate #5 by 10 ys 
and is used as a noise gate 

S&H gates #2 - 20 are contiguous 

S&H gate pulsewidths 200, 400, 800 ns 

Increase of pulsewidth is controlled by 
.comparison of outputs of S&H gates #1 and #20 


63 



TABLE 13 (cont.) 

Decrease of pulsewidth is controlled by 
comparison of outputs of S&H gates #1 

Bandwidth of pulsewidth decision 

Number of 10-dB channels 

AGC bandwidth 

Instantaneous dynamic range of 
received power 

Total dynamic range 

Footprint 

Fresnel Zone 

Diffuse return (max. ) 


and #10 
1 Hz 
3 

10 Hz 

30 dB 
60 dB 

90 m 
7 km 


The initial range of the radar return is acquired and 
tracked by a comparison of the signal with a delayed version 
of the waveform in a threshold device, which also provides a 
timing pulse for positioning S&H gate #5. The 20 S&H gates 
are spaced as indicated in Table 13. Gate #1 is used to 
monitor the noise level of the return signal and to provide 
a reference for the other S&H gates so that the full return 
signal is always captured. The outputs consist of those of 
the three sets of 20 S&H gates, which correspond to three 
power levels, plus that of the initial threshold detector. 

This satellite terrain radar system should be basically 
used to construct a radar topographic map of the earth, 
which would also include other terrain intelligence such as 
vegetative cover and dielectric constants from calibrated 
specular-reflection data. Accurate height measurements over 
flat areas can provide continental networks of reference 
points, which can be used to connect intercontinental 
systems and derive for the first time relative orthometric 
heights on a global scale. A base map of this type, once 
established, can then be used to determine future changes in 
terrain features caused by climate, weather, or man-made 
effects. 

5.4 Conclusions 

A large data base for terrain radar altimetry was 
obtained by the altimeter aboard the Sky lab mission. 

Although there was considerable latitude in the quality of 
the data due to complexities of the terrain and limitations 
in instrument performance, the latter consisting primarily 
of antenna pointing error and bandwidth restriction, -the 
available data base showed the existence of a distinct radar 
signature of the low-frequency profile of the elevations 
along the ground track. Different levels of return signal 
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strength and different signal wave forms were identifiable 
with gross terrain features. Loss of lock in range tracking 
and consequently loss of significant terrain intelligence 
over land areas of complex relief indicated the need for 
larger dynamic power range and range rate capability in a 
dedicated terrain altimeter, and the specifications for such 
a system have been outlined. 
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FIGURE 1. GEOMETRY OF FLAT TERRAIN OBSERVATIONS. 
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FIGURE 2, GEOMETRY OF COMPLEX TERRAIN OBSERVATIONS, 
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FIGURE 5. ACQUISITION MODE, 





FIGURE 6. SAMPLING GATES FOR RANGE-TRACKER AND WAVEFORM , 
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FIGURE 7. RANGE-TRACKER GATE AND CONTRul CURVE. 



FIGURE 8. DATA FLOW DIAGRAM, 













FIGURE 9. NRL DATA REDUCTION AND ANALYSIS 
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FIGURE 10, SKYLAB ALTIMETRY GEOMETRY. 




FIGURE 11. SKYLAB ALTIMETER FOOTPRINT 






FIGURE 12. SKYLAB GROUNDTRACKS 
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FIGURE 13, FINE-STRUCTURE AGC AND ALTIMETRY GROUND HEIGHT (PASS 16). 
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PHOTOGRAPHY, ALTIMETRY GROUND HEIGHTS 









PHOTOGRAPHY, ALTIMETRY GROUND HEIGHTS, AND AGC (SL-3, PASS 28) , 










0“) 1H9I3H 

aNnoyo 

AM13WI11V 


SKYLAB PHOTOGRAPHY, ALTIMETRY GROUND HEIGHTS, AND AGC (SL-3, PASS 28) 







FIGURE 21. COMPARISON OF SKYLAB ALTIMETRY GROUND HEIGHTS AND TOPOGRAPHIC 
MAP ELEVATIONS (SL-2, PASS GT20) . 
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FIGURE 22. COMPARISON OF SKYLAB ALTIMETRY GROUND HEIGHTS AND TOPOGRAPHIC 
MAP ELEVATIONS (SL-3, PASS 16) . 






MAP ELEVATIONS (SL-3, PASS 18). 







FIGURE 25. COMPARISON OF SKYLAB ALTIMETRY GROUND HEIGHTS AND TOPOGRAPHIC 
MAP ELEVATIONS (SL-3, PASS 28). 





FIGURE 26, AGC AND TOPOGRAPHIC HEIGHT RESIDUALS (SL-2, PASS 1), 









FIGURE 27, AGC AND TOPOGRAPHIC HEIGHT RESIDUALS (SL-3, PASS 16). 
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FIGURE 28. AGC AND TOPOGRAPHIC HEIGHT RESIDUALS (SL-3, PASS 17). 









100 





IDAHO WYO 


4ui)_MV-S1VnaiS3a±H913H 


o 

in 


0 
in 

1 


o 

O 


O 

in 


I 


I 


o 

0 

C\J 

1 


o 

o 



(W0P) 09V 




4 

..jJHI 





FIGURE 32. OCEAN RADAR RETURN PASS 11 (SL-3) DB 707-708, 759-760, 815-816 
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FIGURE 39. TERRAIN WAVEFORMS, RAMP AND PLATEAU GATES (SL-3, PASS 18) 
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FIGURE 41. TERR 




FIGURE 42. DISTRIBUTION OF RADAR CROSS-SECTION. 
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FIGURE 43. EXAMPLE OF TOPOGRAPHIC MAP PROFILE (SL-3, PASS 16). 
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FIGURE 45, MEAN AND STANDARD DEVIATION OF HEIGHT DIFFERENCE 
AND CROSS-CORRELATION OF ALTIMETRY GROUND HEIGHT 
VERSUS MAP ELEVATION (SL-3, PASS 16, DB 204-228) 
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FIGURE 46. SKYLAB ALTIMETER OUTPUT AND TOPOGRAPHIC MAP PROFILE (SL-3, PASS 16, DB 246-270). 
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FIGURE 48. WAVEFORM AND STATISTICAL CHARACTERISTICS OF RADAR RETURN. 








FIGURE 49, RAMP GATE (NO. 7) OUTPUT BEFORE RANGE-TRACK LOSS (DB 243). 
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FIGURE 51. ACQUISITION AND TRACKING SYSTEM FOR A DEDICATED TERRAIN SATELLITE 
RADAR ALTIMETER. 






























APPENDIX A 


CATALOGUE OF TERRAIN SIGNATURES 

KEY TO COLUMN HEADINGS 

Column 

1. Data Block 

2. • SM and S M - sub-mode and sub-sub-mode of operation 

(see Section 1.7). 

3. AGC (dBm)-AGC signal level in dB above mW, calibrated. 

4. Rise-Time (Gates) - Number of S&H gates required to 
cover 10%-90% rise of amplitude of ramp portion of 
signal. Note S 2 M-1 operation. When the signal does 
not rise from the noise level, the number in parentheses 
indicates the amplitude of the signal in the first gate 
as a percentage of the maximum amplitude. This is an 
indicator of the slope of the Signal ramp, a long 
rise-time signifying complex terrain and reflectivity 
characteristics. 

5. Auto-Corr. - Normalized auto-correlation function value 
for 100 samples of No. 4 S&H gate signal with jump index 
or lag of 0.01 seconds (72 meters). 

6. h * { ra )- Reduced altimetry ground height, _ corrected for 

m 

bias for each pass. In general, this correction is 
obtained from a water surface or a large homogeneous 
terrain area such as a desert, where absolute height 
can be easily derived from a topographic map. 

Preceding page blank not filmed a_1 j 


7 . 


h (m) - Average of maximum and minimum map elevations 
mid * 

in a footprint of a data - 

^h/2 (m) - Absolute value of one-half of the difference 
--- between the maximum and minimum elevations within a 
footprint of a data block. 

9. Min. Res (m) - Defined as MIN. RES. = h m * -(h mid +A h/ 2 *-° 
or MXH.RBS. - V ' (h mid ~ L h/2 ) - °* 

10. Type of Terrain - Self-explanatory. 

11. Location - Self-explanatory. 
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Catalogue of Terrain Radar Signatures, SL-2, Pass 1, 30 May 1973 
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Table A- 2 (continued) 
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TABLE A- 3 

Catalogue of Terrain Radar Signatures, SL-3, Pass 17, 7 September 1973 
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TABLE A- 3 (continued) 
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TABLE A- 3 (continued) 
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